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■ We propose a model for the hard X-ray (> 10 keV) emission observed from the 
supernova remnant Cas A. Lower hybrid waves are generated in strong (mG) magnetic 

^ . fields, generally believed to reside in this remnant, by shocks reflected from density 

, inhomogeneities. These then accelerate electrons to energies of several tens of keV. 

^] \ Around 4% of the x-ray emitting plasma electrons need to be in this accelerated distri- 

00 I bution, which extends up to electron velocities of order the electron Alfven speed, and 

' is directed along magnetic field lines. Bremsstrahlung from these electrons produces the 

■ observed hard x-ray emission. Such waves and accelerated electrons have been observed 

■ in situ at Comet Halley, and we discuss the viability of the extrapolation from this case 
to the parameters relevant to Cas A. 
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I Subject headings: acceleration of particles — radiation mechanisms: non-thermal — shock 

^ ' waves — supernova remnants 
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^ . 1. Introduction 

A number of supernova remnants have recently been recognized as likely acceleration sites for 
cosmic ray electrons. SN 1006 (Koyama et al. 1995) and G347.3-0.5 (Slane et al. 1999; Koyama et 
al. 1997) both show spectra with continua so strong in the spectral region around 1 keV that no line 
emission is detectable with e.g. the Advanced Satellite for Cosmology and Astrophysics (ASCA). 
Laming (1998) has shown for SN 1006 that any explanation based on non-thermal bremsstrahlung 
fails. It must either produce line emission that should be observable, or in the most extreme case 
of a shock propagating through a pure carbon plasma (which will produce no lines in the bandpass 
of instruments such as ASCA), the amount of carbon required to produce the observed intensity 
is well in excess of that expected for a Type la supernova, and indeed exceeds the Chandrasekhar 
limit by large factors. Thus synchrotron radiation from cosmic ray electrons with energies of order 
100 TeV is left as the only viable interpretation. This conclusion is essentially confirmed by the 
detection of TeV 7-rays from SN 1006 by Tanimori et al. (1998). 
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A second class of cosmic ray accelerator candidates do show the line emission expected, but 
also have a hard X-ray "tail" to their X-ray spectra. Members of this class include Cas A (Allen et 
al. 1997; Favata et al. 1997), IC 443 (Keohane et al. 1997), and possibly RCW 86 (Allen ct al. 1998). 
While less compelling than the first class, no other mechanism has been postulated to date that 
could explain the approximate power law behaviour of the emission, seen in Cas A up to energies in 
excess of 100 keV (Allen et al. 1997; The et al. 1996). It is the purpose of this paper to investigate 
such a possibility, where collective motions of the plasma (specifically lower hybrid waves) excited 
by ions reflected from shock front decay accelerating electrons into a non-Maxwcllian distribution 
function. We begin by reviewing the simplest models for collisionless shocks. Section 2 reviews 
the theory of electron energization by lower hybrid waves, which has received recent attention 
in connection with cometary X-ray emission (Bingham et al. 1997; Shapiro et al. 1999). Section 
3 describes the application of these results to Cas A, while section 4 discusses the assumptions 
underlying this approach, and some of the consequences for our understanding of Cas A that such 
a model, if correct, would yield. 

AssTiming that the gas pressure >> magnetic pressure in the downstream plasma, the temper- 
ature Ta reached by particle species a with mass rua is given by conservation of energy, momentum 
and particles across the shock discontinuity by 

3 

ksTa = Y^ruav'^^ (1) 

for shock velocity Vg ■ Thus in the limit of true collisionless plasma, the particles will have temper- 
ature proportional to their masses, i.e. protons will have a temperature mi/rrif. = 1836 times that 
of the electrons. Obviously in realistic conditions Coulomb equilibration will take place, with rate 

^-_0 13^^ (2) 

in c.g.s. units, where AT = Ti — Tg, and the ions have charge Z and atomic mass A. So long as 
Tj » Tg we can integrate to get 

Te = f - X O.U-TiUet] . (3) 

Averaging this temperature over the life of the shock modifies Tg by an extra factor 5/7. Such a 
model, assuming values of Z in the range 8-26 (corresponding to oxygen and its burning products) 
adequately explains the observed temperatures of shocked ejecta and circumstellar medium (CSM), 
1.25 and 3.8 keV respectively for UgT ~ 10^^ cm~^s in each case in Cas A (Favata et al. 1997), but 
cannot explain the hard x-ray continuum extending out to 100 keV. 



2. The Electron Distribution Function from Lower Hybrid Wave Damping 



From observations of radio synchrotron emission Cas A is believed to have relatively high 
magnetic fields, perhaps up to the equipartition value of order 1 mC (Gull 1975; Keohane, Gotthelf, 
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& Petre 1998). A total energy in cosmic ray electrons and magnetic field of ~ 10 erg is required 
(c.f. Longair 1994), which for maximum synchrotron luminosity must be equipartitioned between 
the two. If distributed uniformly throughout the Cas A shell, this evaluates to a magnetic field 
of 0.5 mG. Thus lower hybrid waves excited by cross field ion motions can accelerate electrons 
into a high energy tail onto the otherwise thermal Maxwellian distribution function. These are 
electrostatic ion waves that occur when the electrons are magnetized, which is likely to be the case 
at so-called quasi-perpendicular shocks. Due to the the shock compression of the magnetic field 
component perpendicular to the shock velocity vector, quasi-perpendicular shocks are likely to be 
more prevalent than quasi-parallel shocks in Cas A, and in supernova remnants in general, and will 
be assumed throughout this paper. Electron velocities at least up to the electron Alfven velocity, 
B I yf^Kn^rrie where B is the magnetic field strength, and ng and me are the electron density and 
mass respectively, are possible. This corresponds to energies of order 100 keV for Cas A, making 
it an attractive mechanism to explore. 

Electrons accelerated by these waves emit bremsstrahlung radiation, and have been recently 
discussed as a mechanism for comctary x-ray emission. Neutral atoms or molecules outgassing from 
a comet nucleus can be photoionizcd by solar UV radiation. In the interaction of these stationary 
(in the comet rest frame) ions with the ambient solar wind, a modified two-stream instability 
develops that in the case of magnetic field perpendicular to the relative comet-solar wind velocity 
generates lower hybrid waves (Bingham et al. 1997; Shapiro et al. 1999). 

We will apply a similar model to Cas A, following in part the analytic treatments of Vaisberg 
et al. (1983), Krasnosel'skikh et al. (1985), and Begelman h Chiueh (1988), where ions reflected 
from shock fronts in Cas A move back upstream through the preshock plasma generating a similar 
modified two-stream instability. Such reflected ions and associated wave generation are common 
features in in situ observations of solar system shocks (c.f. Lengyel-Frey et al. 1997), and of models 
of electron energization at collisionless shocks (c.f. Cargill & Papadopoulos 1988; Tokar et al. 
1986; McClements et al. 1997; Vaisberg et al. 1983; Krasnosel'skikh et al. 1985). Vaisberg et al. 
(1983) give an analytic form for the electron distribution function following lower hybrid wave 
energization, which agrees well with the particle-in-cell simulations performed by Shapiro et al. 
(1999). The mathematical development of the theory for these waves is summarized in Appendix 
A. It is shown that for both reactive and kinetic forms of the instability, wave growth is strongest for 
wavevectors with a small parallel component h^^jh < uOpi/uOp^. Hence the wave can simultaneously 
be in resonance with ions moving across and electrons moving parallel to the magnetic field, since 
the phase velocities are much larger along the magnetic field than across it. This facilitates energy 
transfer between electrons and ions. The resulting electron distribution function parallel to the 
magnetic field is also derived in Appendix A. 

However the wave group velocity away from the shock front is lower than the shock velocity 
itself for k\\/k < ujpi/cope- Only at k^\/k = Upi/upe, or greater can the waves stay ahead of the shock. 
This is important because the growth rates, particularly for the kinetic form of the instability, are 
not large, and for smaller A;|| the waves are likely to be overrun by the shock before significant 
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growth can occur. Hence in the following we assume that in the preshock region, a significant level 
of wave activity exists only for k\\/k = uj^ijujp^. With this simplification the accelerated electron 
distribution function is then given by equation (A17) with cos^ = ojpi/ujpe 




where VAe = VAi sj rniUi / niene is the electron Alfven velocity, (vai = B/y/ATrriimi being the usual 
Alfven velocity), is the density of shock reflected ions, and is a constant of integration, 
denoting the maximum electron velocity. In the quasi-linear theory used in Appendix A, the 
accelerated electron density, n'^, or Vm remains a single free parameter. They are related by 

where an extra factor of two has been included to include the electron velocity range —Vm — ^ +Vm- 

For a reflected ion population ~ 25% of the preshock ion density (c.f. Cargill &; Papadopoulos 
1988; McClements et al. 1997, and references therein), we get 

^ = 0.06^4 + ^1-4- (6) 

Anticipating from the following section that n'^/rie — 4%, vac ~ 26vs and vn = 0.53^^ (assuming 
complete equipartition of magnetic and ion thermal energy, but insignificant electron-ion equilibra- 
tion) gives v„i = VAe Vvi + ^ — 1 where A is the atomic mass of the ions. Hence Vm ~ vac, slowly 
increasing with the ion mass. 



3. Application to Cassiopeia A 

We compute model bremsstrahlung spectra using the distribution function given by equation 
(4) above for the electron velocity parallel to the magnetic field, representing the electron velocity 
perpendicular to the magnetic field as Maxwellians. The temperature for these perpendicular 
components is taken to be the same as that fitted to the thermal bremsstrahlung emitted at 
energies less than 10 kcV. A fit to the BcppoSAX Medium Energy Concentrator (MECS Boella 
et al. 1997) spectrum from 1997 November 26 determined this temperature at 3.3 keV. This fit is 
shown in Figure 1. Other workers find temperatures in the range 2.9-4.2 keV (Allen et al. 1997; 
Favata et al. 1997; Vink, Kaastra, &: Bleeker 1996) from ASCA and BeppoSAX data. 

The non-relativistic bremsstrahlung cross section da for photon emission with energy in the 
range E ^ E + dE is given by (Berestetskii, Lifshitz, k. Pitaevskii 1982) 

^ 647r2^2 2 "^eC^ Pf 1 dE 

"""^ {Pi - Pf? ^ (1 - exp {-^-kZIv})) (exp {2'iTZ/vi) - 1) £^ 

(7) 



da = — 



dz 



^ , iZ iZ 
2F1 —,—,1,2 

Vf Vi 



-5- 



where pi, pf, Vi, Vf are initial and final electron momenta and velocities, Z is the nuclear charge, a 
is the fine structure constant, and rg is the classical electron radius. Using standard formulae for 
manipulating the hypergeomctric function 2F1 and its derivative in the limit of low Z/vi and Z/vf 
a simpler form suitable for evaluation of the double integral to compute the emission rates can be 
found; 

da = ^Z'ar'-^^I^ ^/^» W ( EllPl] ^ (8) 

3 ''vf{l-exp{-27rZ/vf)){exp{27rZ/v,)-l)^^\p,-pfJ E' 

In checks against bremsstrahlung spectra computed with the full formula given by equation (7) 
(Sutherland 1998), the simpler equation (8) was found agree at the few percent level for Maxwellian 
electron distributions, for temperatures appropriate to Cas A for target charges up to 8 (i.e. a fully 
ionized oxygen plasma). Cas A has very little hydrogen. Oxygen is probably the most abundant 
element (Vink, Kaastra, & Bleeker 1996). 

In Figures 2 and 3 we plot spectra from the BcppoSAX MECS and the higher photon energy 
Phoswich Detector System (PDS Frontcra ct al. 1997) instruments, together with model spectra. 
The softest model spectrum in Figure 2 is a pure thermal bremsstrahlung spectra at 3.3 kcV in 
fully ionized oxygen. The emission measure for this spectrum is UeUoV = 3.05 x 10^^ cm^'^ for 
an assumed distance to Cas A of 3.4 kpc (Reed et al. 1995). The successively harder spectra have 
been computed by taking 4% of these thermal electrons and putting them into a lower hybrid 
wave energized distribution function, for electron Alfven velocities of 32, 40, and 48 atomic units 
(in atomic units, e = h = rrie = 1, so the atomic unit of velocity is 2.188 x 10^ cm s^"*^), with 
Vm = 1.77f^e- Figure 3 shows the same data with the pure thermal bremsstrahlung and non- 
thermal bremsstrahlung models for vac of 56, 68, and 80 atomic units for a pure He plasma, with 
Vm = l.llvAe- In both cases, the curves with Vm ~ 70 atomic units appear to give the best match 
to the data up to energies ~ 50 keV. 



4. Discussion 

4.1. The Morphology of Cas A 

One immediate problem is that the niG magnetic fields are present only inside the Cas A rem- 
nant (presumably amplified by Rayleigh- Taylor instabilities at the contact discontinuity), whereas 
most shock models require such fields to be present ahead of the shock (to be discussed in more 
detail below). However Cas A is an extremely inhomogeneous supernova remnant. It has numerous 
knots of optical emission from low charge states, which must be considerably denser than the sur- 
rounding plasma so that the recombination and radiative cooling times are sufficiently short. The 
blast wave travelling through the circumstellar material will to split into transmitted and reflected 
shocks upon meeting these inhomogeneities (see Sgro 1975; Borkowski, Blondin, Sz Sarazin 1992, for 
fuller discussion). Indeed exactly this happens in the simulations of Borkowski et al. (1996), giving 
rise to the possibility that shock waves are propagating throughout the X-ray and radio emitting 
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shell. Since the plasma has already been heated by the blast wave or reverse shock these shocks 
will most likely be of relatively low Mach number (but still supercritical; see subsection 4.5 below). 
Another possibility for inhomogcncous or clumpy SNR is that the Raylcigh- Taylor fingers from the 
unstable contact discontinuity can propagate outwards and penetrate through the forward shock 
(Jun, Jones, & Norman 1996). The by-now famous Chandra x-ray images (Hwang, Holt, k, Petre 
2000) suggest that both scenarios could indeed be true. The forward shock is visible as a faint ring 
around the outside of the bright x-ray shell, which has its origin in shocked ejecta. The boundary 
between faint and bright emission is the contact discontinuity, which has no obvious penetration 
through the blast wave, except in the "jet" region in the NE quadrant. However the apparent "in- 
version" of the presumed pre-supernova element stratification (Hughes et al. 2000; Hwang, Holt, & 
Petre 2000), with Fe being found in outer regions of the eastern part of the remnant with respect to 
lighter elements like Si, suggests a role for the Raleigh- Taylor instability, and makes it more likely 
that these regions have come into contact with strong magnetic fields at the contact discontinuity. 
Lower hybrid electron acceleration is also more likely in the heaviest element plasmas, (i.e. the Fe 
ejecta), for reasons discussed in the next subsection. 



4.2. The Thermal Electron Temperature 

The lifetime of accelerated electrons against Coulomb collisions with the thermal background 
is t ~ 5 X 10^ {E/1 keV)'^''^ s which evaluates to 14 - 150 years for electron energies of 20 - 100 keV. 
Hence electron acceleration must be ongoing, or have ceased in the very recent past, though this 
last possibility seems unlikely since the hard x-ray emission from Cas A has been known since the 
earliest days of x-ray astronomy (Gorenstein, Kellogg, k. Gursky 1970). 

In section A.4 of the Appendix it is shown that an appropriate criterion for the presence of 
lower hybrid waves is that the electron gyroradius should be less than the wavelength, which can 
be expressed as Tg/Tj < Vt^^jZ juj^^. In a magnetic field of 1 mC and electron density of 10 cm~^ 
this evaluates to Tg << 3 x 10^ (^;c,./4000km s^"*^)^. The shock velocity has been determined from 
proper motion studies. For the same assumed distance as before, radio data give an expansion 
velocity of ~ 2000 km s~^ (Anderson h Rudnick 1995), while X-ray studies give velocities of ~ 
3500 km s~^ for the bright ring at a radius of 110", and ~ 5200 km s^^ inferred for the blast wave 
at a radius of 160" (Vink et al. 1998), all of which require T^, significantly lower than that observed. 
As discussed in section A.4, an alternative condition on the electron Landau damping rate requires 
/c^T << meJll^/A;^, which is the same as above for k^^jk ~ ujp-Jujpf,, but less restrictive for smaller 
k\\. The first criterion is also not met in the plasma around the nucleus of Comet Halley where lower 
hybrid waves and accelerated electrons were directly detected in situ (Gringauz et al. 1986; Klimov 
et al. 1986). They observed VLlh ~ 100 rad s^S and Ue 300 cm^^ giving O^/cj^g ~ 2 x 10"^, 
while Te ~ 2 X 10^ K and Tj ~ 10^ K. 

One plausible solution to this dilemma is that in the region of wave generation, the electron 
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temperature is in fact lower than the estimates and observations above would indicate. This is 
possible if the electrons are not energized on passage through the shock due to being bound to ions, 
only becoming ionized some time later on either by collisions in Cas A or by photoionization in 
Comet Halley. This is demonstrated in calculations of the Cas A reverse shock propagating through 
clumps of Fe ejecta by Mochizuki et al. (1999), where temperatures satisfying the condition on the 
electron gyroradius are found. Thus for Cas A the hard x-ray emission is more likely to be associated 
with regions of shocked ejecta where many electrons are still being ionized from bound states. In 
order that this putative ionizing region in Cas A should produce enough bremsstrahlung emission 
to be visible against the thermal continuum, a higher background ion charge might be necessary. 
This is also more likely if the thermal bremsstrahlung continuum with Tg = 4 x 10^ K comes from 
the shocked circumstellar medium, composed of He, N, and possibly O, while the non-thermal 
component comes from the ejecta, composed of O and heavier elements. Thus the local fraction of 
accelerated electrons may be different from the 4% determined above by matching to the BeppoSAX 
MECS/PDS spectra. However the maximum accelerated electron velocity Vm oc (n'^/rie)^^^ is 
relatively insensitive to this, making large changes in the bremsstrahlung spectrum from that 
calculated before unlikely. The temperature assumed above for the perpendicular components of 
the accelerated electron distribution will also change. Again this will not significantly alter the 
non-thermal bremsstrahlung spectrum. 

Another reason for spatially separating the thermal and non-thermal bremsstrahlung emitting 
regions is that over the 300 years or so of the remnant's existence, collisions between the accelerated 
and thermal electrons would be likely to have raised the thermal electron temperature to a value 
well in excess of that observed if they were cospatial. 



4.3. The Electron Alfven Speed 

Another crucial parameter in determining the visibility of non-thermal bremsstrahlung is the 
electron Alfven velocity. From observations of radio synchrotron emission, Cas A is believed to 
have a large magnetic field ~ 1 mC, generated by turbulence associated with Ray leigh- Taylor 
instabilities at the contact discontinuity between the shocked ejecta and circumstellar material (c.f 
Gull 1975; Keohane, Gotthelf, & Petre 1998). Assuming such instabilities generate magnetic fields 
up to equipartition, we may write 

1 

— = - {niksTi + neksTe) (9) 

where Tj and Tg are given by the shock jump conditions (equation 1). The resulting electron Alfven 
velocity thus becomes 

VAe = \ -r~ = Vs\ = 2Qvs (10) 

y 47rneme V o'me 

for fully ionized plasmas where nuclear mass Arup is such that A = 2Z. The X-ray blast wave veloc- 
ity (discussed above) gives VAe — 60 atomic units according to equation (10), with correspondingly 
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lower values from the other observations. This is more than adequate to give a sufficiently high VAe 
for the non-thermal bremsstrahlung to originate in oxygen or heavier element plasma. Hard x-rays 
from a helium plasma require higher vac, still consistent with the observed x-ray blast wave, but 
not with the lower inferred shock velocities. 

4.4. High Energy Behaviour of the Continuum and ^^Ti 

The highest energy data points (above 100 keV) may still possibly require the existence of 
cosmic rays. As can be seen in Figures 2 and 3, vac or cannot be arbitrarily increased to repro- 
duce this emission without overestimating the non-thermal continuum at lower energies. However 
in principle relativistic electron energies are possible from lower hybrid wave acceleration (Mc- 
Clements et al. 1997; Vaisberg et al. 1983). The maximum cutoff in our approach stems from 
taking only waves with k^\/k = tOpi/uJpe, which satisfies the requirement that the perpendicular 
group velocity of the waves be greater than the shock velocity, to allow sufficient time for wave 
growth. With stronger reactive instabilities, waves with ^ may exist at sufficient intensity to 
make electron acceleration by lower-hybrid waves a plausible mechanism for solving the so-called 
electron injection problem (McClcmcnts et al. 1997). This arises because first order Fermi accel- 
eration at quasi-perpendicular shock fronts is only effective for electrons of mildly relativistic or 
higher energies (Levinson 1996), and so some other "injection" mechanism is required. 

The question of the spectral behaviour at high energies is further confused slightly by the 
discrepancy in fluxes observed above 100 keV by BeppoSAX and OSSE. Compared with the data 
from OSSE presented in Figure 2 of Allen et al. (1997), the BeppoSAX PDS fluxes are up to an 
order of magnitude higher above 100 keV. Cas A is unlikely to be a variable source at these (or 
any other) energies, and so this problem presumably lies with the calibration of one or other of the 
instruments. 

Another issue complicating the comparison of continuum spectra, but of supreme interest in 
its own right, is the possibility of detecting 7-ray line emission from the nucleus ^^Sc at energies of 
67.9 and 78.4 keV. ^^Sc is a decay product of ^^Ti, which is one of the most important radioactive 
elements produced in the a-rich freeze out of core collapse supernovae, its abundance also being 
sensitive to the position of the mass cut (c.f Vink et al. 2000, for further discussion). The break 
in the lower-hybrid wave accelerated electron spectrum occurs at energies very close to where the 
^^Sc lines should be present, and the substitution of the power law in the fit to the BeppoSAX PDS 
data by Vink ct al. (2000) with the non-thermal bremsstrahlung spectrum modeled in this paper 
might lead to interesting results for the flux in these lines. 
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4.5. Properties of CoUisionless Shocks 

So far we have invoked the excitation of lower hybrid waves ahead of a shock front, and modeled 
the electron distribution function on the assumption that a quasi-steady state can be reached in 
the shock rest frame. We will discuss these assumptions in more detail here. 

In a number of models, various types of plasma waves are excited by ions that are reflected 
from the shock front and travel back upstream. This generally arises in perpendicular shocks due 
to the cross shock ambipolar electric field in the magnetic field overshoot region. The simulation 
of Cargill & Papadopoulos (1988) has reflected ions exciting electron plasma waves and then ion 
acoustic waves to provide eventual electron heating to around 20% of the postshock ion temperature. 
McClements et al. (1997) looked at the excitation of lower hybrid waves by ion gyrating in the 
magnetic field upstream of a fast shock, with a view to modeling the injection mechanism for the 
Fermi acceleration of electrons to cosmic ray energies at perpendicular shocks. A certain amount of 
support for these models comes from observations of shock waves in the solar wind where electron 
plasma and ion acoustic waves are seen ahead of shock fronts (Lengyel-Frey et al. 1997), and lower 
hybrid waves (Klimov et al. 1986) and their associated electron distribution function (Gringauz et 
al. 1986) have also been observed at comet Halley, with Shapiro et al. (1999) recently demonstrating 
by particle-in-cell simulations that the observed wave activity is indeed consistent with the observed 
electron distribution function. 

An important point concerns the nature of the two-stream instability. All of the models 
discussed above assume the ions refiected from the shock to be essentially monoenergetic, thus giving 
rise to fast growing reactive instabilities. At perpendicular shocks this is necessary since the reflected 

ions only move back upstream for about an ion gyroradius before begin swept back behind the shock 
front, thus allowing only a finite amount of time for the instability to develop. This monoenergetic 
property of the reflected ions surprisingly appears to persist in shock simulations to high Mach 
numbers, where the shocks themselves should be highly turbulent (Cargill Sz Papadopoulos 1988; 
Tokar et al. 1986). Both sets of authors comment that the ion reflection in their simulations 
becomes "bursty", giving a time average of 25%; not too different from the fraction of preshock 
ions observed to be reflected at lower Mach number (less turbulent) shocks in the solar wind or 
laboratory experiments. Woods (1987) emphasises that such ion reflection is characteristic of a 
very thin shock transition. For thicker shocks, the ions will gyrate in the magnetic fleld before 
reflection to varying degrees, giving rise to a reflected ion distribution still monoenergetic, but 
having a flnite angular spread. Woods (1987) terms this a transition from "specular" to "diffuse" 
reflection. Gedalin (1996) has considered this distinction in more quantitative detail. Diffusely 
reflected ions will produce slower growing instabilities ahead of the shock, and hence produce less 
electron heating, and have been suggested as a possible reason by Laming (1998) for the low amount 
of collisionless electron heating observed at fast (~ 2600 km s^^) shocks in the NW region of SN 
1006 by Laming et al. (1996). 

We have argued above that based on our knowledge of the properties of Cas A, the shocks 
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producing the reflected ions and lower hybrid waves are hkely to be of relatively low Mach number, 

since they arc reflected from density inhomogeneities and are propagating through plasma that 
has already been heated by the forward or reverse shocks. Hence a "bursty" ion reflection is 
unlikely, and the assumption of steady state conditions embodied in equation (A15) remains valid. 
The sound speed in plasma shocked by the blast wave with speed Vs and with no electron-ion 
equilibration is y^5kBTi/3mi = \/h/lQvs- This plasma continues expanding at speed O-TSus, and a 
shock reflected from a density inhomogeneity will hit it with relative velocity -s/lbAr/lQ/ (4 — Ar)iis 
(i.e. a Mach number of ^/'iAj-j (4 — ^r)), where Aj. is related to the density contrast, A, by 



A = 3Ar {4:Ar - 1) / ( y^3Ar {4:-Ar)-V5 {Ar - 1) ) (Sgro 1975). For A ^ oo, Ar = 2.5 and 



the maximum Mach number of the reflected shock is -v/5- In order that these secondary shocks 
propagating back through the shell arc supercritical (i.e. that they reflect ions to form a precursor 
that excites the lower-hybrid waves), they must have Mach number greater than the so-called first 
critical Mach number. Edmiston & Kennel (1984) in a survey of first critical Mach numbers for 
various MHD shocks find a value of 1.7 for perpendicular shocks in plasmas with 7 (ratio of specific 
heats) of 5/3 and /? (gas pressure/magnetic pressure) of ~ 1. This requires A^. ~ 2 and a density 
contrast of ~ 30, with lower values required for higher /3. For reference the electron density in 
the shocked x-ray emitting plasma is in the range 10-15 cm~^ (Vink, Kaastra, &: Bleeker 1996; 
Favata et al. 1997), while that in the shocked optically emitting knots determined from S II line 
ratios is ~ 10^ — 10^ cm"^ (Chevalier &; Kirshner 1978; Reed et al. 1995). The initial mass density 
contrast will be increased from the simple ratio of these electron densities since the x-ray emitting 
plasma is more highly ionized than the optical knots, by a factor of ~ 10. However the density of 
the optical knots may also have been enhanced over that due to shock compression by radiative 
cooling. Chevalier &: Kirshner (1978) report densities of < 100 cm~^ from diffuse (presumably 
unshocked) S II, which given the difference in ionization state between this and the x-ray emitting 
plasma still gives sufficient density contrast to refiect supercritical shocks back into the shell. 

Whether the ion refiection is diffuse or specular is a harder question to answer. We have 
assumed a kinetic form for the instability, since this is likely to be more realistic, and also gives 
a simpler derivation of an analytic form for the electron distribution function. In order to get 
sufficient time for wave growth to occur ahead of the shock, we are limited to waves with parallel 
components of wavevector = ujpi/tOpe or greater, in order to have group velocities away from 
the shock comparable to the shock velocity itself. Specularly reflected ions give rise to much faster 
growing instabilities, where waves with smaller fcjj may grow to significant intensities before being 
overrun by the shock to allow harder electron spectra to be produced, and might conceivably 
produce a better match to the BeppoSAX MECS/PDS data. This issue of specular versus diffuse 
ion reflection and the reality of the electron distribution function discussed here is likely to remain 
an open question at least until high resolution X-ray spectra become available. We note that the 
electron distribution function produced by the wave activity should have discernable effects on 
the line spectrum. Ionization and excitation rates will be altered by the extra population of high 
energy electrons, and modeling of the X-ray spectrum might be the best way to confirm or refute 
the existence of such an electron distribution function (see Gabriel et al. 1991, for an example). 
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Such an effort is beyond tfie scope of this paper, and requires high spectral resolution observations 
where individual lines can observed. 

5. Conclusions 

We have presented a model based on plasma wave excitation by collisionless shocks to gener- 
ate a non-thermal electron distribution, the bremsstrahlung from which can explain the hard X-ray 
"tail" in the spectrum of Cas A, for very reasonable model parameters. This is important because 
the alternative interpretation, that of X-ray synchrotron emission from cosmic ray electrons with 
energies of order 10^^ eV, appears to have been accepted largely on the basis that it is the only 
interpretation available, though some further observational support for this view is accumulating 
(Vink et al. 1999; Hwang, Holt, & Petre 2000). Our model is based on observation of the relevant 
waves and energized electrons at Halley's comet, and analytical and numerical treatments of the 
process confirming the consistency of the observed electron distribution with the observed wave 
activity. We have assumed that the wave excitation is associated with shocks which propagate 
throughout the shell of Cas A, with the waves mainly being generated in ionizing regions of the 
ejecta, although other mechanisms could be viable. For instance, Begelman & Chinch (1988) dis- 
cuss similar waves excited by relative electron-ion drifts in Advcction Dominated Accretion Flows. 
It is also worthwhile to point out that electron acceleration by electrostatic waves is essentially 
instantaneous on the timescales relevant to supernova remnant evolution, whereas for example 
stochastic mechanisms (i.e. second order Fermi acceleration) are not (c.f. Blasi 2000). 

Similar explanations may also be valid for the supernova remnants RCW 86 and IC 443, for 
which hard X-ray tails have also been observed. Chevalier (1999) has also expressed scepticism 
that such emission in IC 443 is due to shock accelerated cosmic ray synchrotron emission. Our 
discussion has concentrated on Cas A simply because of its well known age, distance, shock ve- 
locities, and conspicuous high energy continuum. If proved correct, then this explanation could 
provide the clearest evidence to date for equipartition of magnetic and thermal energy densities, 
and the precise shape of the non-thermal bremsstrahlung spectrum should be rich in diagnostic 
potential for inferring magnetic fields, shock structure and the likely nature of preshock ion driven 
instabilities. 

Finally we note that our non-thermal bremsstrahlung model is qualitatively different to those 
of previous authors (Asvarov et al. 1990; Sturner, Dermer, &; Mattox 1996). Both references 
consider non-thermal bremsstrahlung in X-rays from the same electron population that produces 
synchrotron radiation emission, and find that for the hard X-ray tail to be explained in this way, 
it is likely that too much radio emission would result. Baring, Jones, & Ellison (2000) give a 
fuller discussion of such scenarios. Our model here postulates a separate electron population which 
is non-relativistic; basically just a small distortion of the thermal distribution, which need not 
necessarily be associated with the cosmic ray electron distribution, and so no great change in radio 
emission should result. 
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A. Appendix: Lower Hybrid Waves 
A.l. Dispersion Relation and Reactive Growth Rate 



For cold magnetized electrons the dielectric tensor is 
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where Upe and are the electron plasma and gyrofrequencies respectively. For longitudinal waves 
we put E||k in the wave equation, take the scalar product with k with the magnetic field along the 
z-axis and set the result equal to zero to derive the dispersion relation; 
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which gives 
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To this we can add a term for thermal non-magnetized ions (l — ^ (^ijj/V2kvi)) oo^/k'^vf where (/) (x) 
is the usual plasma dispersion function and Vi = y/kBTi/nii, the ion thermal velocity. In the limit 
a; >> kvi, (x) — > 1 + l/2x^ + . . . and the well known form for the dielectric tensor for lower hybrid 
waves results 
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where ujpi is the ion plasma frequency. This yields the dispersion relation u) ~ uj^^ as cos 
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0—^1. Q^i? is known as the Lower Hybrid frequency and is equal to with 
Qi being the ion gyrofrequency. We now consider a further ion population moving with velocity U 
with respect to the stationary plasma. The dispersion relation becomes 
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For cos 9 1 equation (A5) takes the same form of that for the Buneman instability where 
Langmuir waves are generated by the ion beam (c.f. Sturrock 1994, pp 113). For cos^ a 
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quartic equation similar to that for the Buneman instabihty results, but with the replacements 
ujpe ^LH and ujpi ClLHi^'pi/i^pi, giving a growth rate 

7 = ^^-^f-V^'- (A6) 



no- 



Moving away from cos ^ = 0, we put a; = Olh in the parallel electron term and the ion term to get 

(a; -kUf = . (A7) 

i + -j^sm y j^cos y ^ 

With kU ~ 0,LH and ujpe » ^Ig, the growth rate 7 is 

/\ 1/2 



u-nLH^±il^±i{-] (A8) 
\ni J cos 

Hence for cosO'^ojpi/upe, the growth rate is of order ^lh, falling off as Jlj/ cos^ for larger cos0. 

As discussed by Begelman & Chiueh (1988) an electromagnetic correction is required in the par- 
allel electron term of K^. Since waves propagate within ^Upi/upe of the perpendicular direction to 
the magnetic field, the wave phase speed parallel to the magnetic field can be very large, approach- 
ing or exceeding the speed of light. This correction can be written LOp^ ^pe/ (l +'^pe/^^'^^) 
(Begelman & Chiuch 1988; Melrose 1986, p 28). We also retain further terms in the expan- 
sion; 1/ (w^ - r^e) ~ -1/^e (1 +'^^/^e + •■■)' teim in parentheses ~ 0.1^/^1 ~ 
WpgU^j/JlgC^ ~ ujpg/k'^c^ ■ v^/v^. Here v'j_ = + is the electron velocity perpendicular to the 
magnetic field, with vai the Alfven velocity, Vg the electron thermal velocity, and c the speed of 
light. Taking k ~ Qe/v±, with A:|| and k± being wavevectors parallel and perpendicular to the 
magnetic field, applying these two corrections yields (see e.g. Vaisberg et al. 1983; McClements et 
al. 1997) 

R-L _ 1 , kl ( __}_ \ '^li ^le /I , , ,2 |l,2'2^-l f 

^ - ^ + ^ 1^^ ^Wl + v^Jv%) (1 + c ) . (A9) 

McClements et al. (1997) discuss reactive growth rates derived using similar dispersion relations 
for various reflected ion distributions. 



A. 2. Kinetic Growth Rate and Accelerated Electron Distribution 

In the presence of suprathermal electrons and ions, two extra terms are added to the expression 
for K^, 

+ % I hA^^ . M,3,-; + ^ ! k^^^ p^3,-.^ (AlO) 
n'ek J w U -k-Ve u'-k^ J U) u -k-Vi OVi 
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where the primes denote quantities with respect to the suprathermal particle distributions, rig and 

n' being the suprathermal electron and ion densities, respectively. The overall growth rate for 
lower- hybrid waves is found by setting uj ^lh + ^7 and taking the imaginary parts, using the 
Landau prescription for the terms in equation (AlO). This produces a growth rate 7; 



{n^H - k . .1) d'vl + f (1 + ^iJk'c'-') ||l.„=n,Wfc„ } ■ (All) 



The first term gives wave growth/damping by the suprathermal ions and the second denotes that 
by suprathermal electrons moving along magnetic field lines. The growth rate due to a Maxwellian 
ion distribution with temperature T moving with bulk velocity U may be evaluated (dropping the 
term in w^/fi^ in the perpendicular electron term); 
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The growth rate is maximized for ■ U — ^lh^ = k^ksT /rrii, hence 
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which gives strongest wave growth at the smallest possible wavevector, kmin = '^^ij \fk^Tjmi 
Hence 

7i,max -0.38-2- — -^ ^, (AM) 

''^'«"Ml + i4^cos2 

and can be seen to be largest for cos Q < Upi/upe, in qualitative agreement with the reactive growth 
rate. However the kinetic growth rate is smaller by a factor ~ Ql/cOp^. If the reflected ions have 
a larger temperature T' than the ambient ion temperature T, the growth rate is smaller still by a 
factor ~ T/T'. Returning to equation(All) the accelerated electron distribution function is given 
(with the assumption of steady state conditions, i.e. 7 = 0) by 

With some manipulation cOpJk'^c'^ = ^LH/k'^''^Ai — ^ii/ ^o^^ ^/^y'^Ai = '"y cos^ 6*/^^^, allowing us 
to integrate to get /g; 



f U ) = ^ exp ( -1 ) ^ - V. + ^ I (A16) 
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where VAe = ^Ai \/ vnini / meUe is the electron Alfven velocity, Vm is a constant of integration, the 

maximum velocity where / (vm) = 0, and 9 is the angle between the magnetic field direction and 
the wavcvcctor where the growth rate is a maximum. We may integrate over /g in the range 
—Vm +Vm to obtain an expression for the number density of energized electrons, n'^. The result 



IS 



n'- ^ l\ ^ni / 9 vt, cos^ 6 



n'e = ' ' exp -- ^ + IllLrrUl ) . (A17) 

Vaisberg et al. (1983) appear to make the approximation (l + u^^/k^c?) ^ 2± X — uj^JlP'f? to derive 
a slightly different form for the accelerated electron distribution, where VAe becomes a firm upper 
limit to v\\. 



A. 3. Group Velocity 



Prom equation (A4) with iOp^. » 
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Hence the wave group velocity away from the shock front is 
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If uj/k± ~ U ~ 2vs then Vg > Vg for k\\/k± ~ cos 6 ~ Upi/upg. Waves excited with lower cos^ 
have lower group velocities. Since the reflected ions only move a distance downstream on the 
order of a gyroradius, growth rates need to be several times Jlj if significant wave growth for these 
waves is to occur before they are overrun by the shock. Waves at cos^ ^ ujpi/upe and above can 
stay ahead of the shock and can hence grow to large intensities, even if the growth rate is small. 
Hence in equation (A17) we simply put cosO = LOpi/ujpe and proceed with the calculations of the 
bremsstrahlung continuum. Similar conclusions result from the use of the dispersion relation with 
electromagnetic and thermal corrections (McClements et al. 1997), and will be discussed in more 
detail in future work. 



A. 4. Applicability Criteria and Connection to Electron Acoustic Waves 

Equation (A16) for the accelerated electron distribution assumes waves are generated at a pitch 

angle to the magnetic field direction given by cos 9 = ujpi/ujpe only, yielding a maximum electron 
velocity of VAe- In realistic conditions both the reactive and kinetic instabilities will excite waves 
with pitch angles down to cos 6* = 0. Consequently waves with phase speeds along the magnetic 



-16- 



field, ^LH/k\\, up to c will exist, allowing some electrons to be accelerated up to relativistic energies. 
However a break in the electron spectrum should still exist at the electron Alfven speed, VAe- 

In the lower hybrid wave the damping of the ion oscillation by electron screening is inhibited 
since the electrons are magnetized. This suggests a criterion for the existence of lower hybrid 
waves of electron gyroradius << wavelength, or v^/^e « Vi/ujpi. This can be reexpressed as 
Te « ^TiO.l/L0p^. Assuming Tj = ^rriivj/kB for a 4000 km s~^ shock in Cas A (assuming 
equipartition between ions and magnetic field) gives Tg « 3 x 10^ K. On these grounds we might 
expect lower hybrid wave acceleration to have ceases in Cas A, since the electrons have been heated 
by coUisional processes to Tg ~ 4 x lO^K. On the other hand this criterion appears not to be met 
for the plasma conditions near the nucleus of comet Halley, where lower hybrid waves and their 
associated electron distribution function were in fact observed in situ in 1986 (Gringauz et al. 1986; 
Klimov et al. 1986), suggesting possibly that it is somewhat naive. From the condition that the 
electron Landau damping rate be << Qlh we derive ksTe << me^\H / k"^]^, which is equivalent 
to the first inequality if k\\/k = cosO = ujpi/ujpe^ but extends to higher electron temperatures for 
cos 9 -^■ 0, which would be more consistent with observations. Shapiro et al. (1999), motivated 
in part by the discovery of cometary x-ray emission, (c.f. Bingham et al. 1997) have recently 
demonstrated by particle-in-cell simulations that the observed wave activity at comet Halley is 
indeed consistent with the observed electron distribution function. 

If in equation (A4) we put 4){x) — > 2x^ — + ... appropriate for hot ions, then with 

fie >> '^pe 

^'-l + 4-$|-0 (A20, 

which is equivalent to equation (3.11) of Begelman &; Chiueh (1988), their dispersion relation for 
electron acoustic waves, in the limit of zero ion drift velocity. These authors derived their dispersion 
relations from fluid equations, assuming the electrons to be strongly magnetized, i.e. Jig > ^pe- 
This requirement is not necessary here. As discussed by Bcgchiian & Chiueh (1988) these electron 
acoustic waves have very similar properties to lower hybrid waves. The growth rate for these waves 
due to a drifting Maxwellian ion distribution may be evaluated in the same way (putting uopi — > 
and — > Upi in the forgoing) with the result 

li,max = Ixl^exp (-1) , /2'2 n \ ^^=F=] ■ (^21) 
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Fig. 1. — Data and fit for the spectral region around the Fe K feature at 6.7 keV taken from the 
BeppoSAX MFCS spectrum. The temperature fitted to the thermal bremsstrahlung continuum is 
3.3 keV. 

Fig. 2. — Data from the BeppoSAX MFCS and PDS instruments, with model bremsstrahlung 

spectra for a fully ionized pure oxygen plasma. The softest model is pure thermal bremsstrahlung 
at 3.3 keV temperature. The harder spectra come from electron distribution functions appropriate 
to lower hybrid wave energization, for electron Alfven velocities of 32, 40, and 48 atomic units. The 
fraction of electrons in the non-thermal distribution is 4%. 

Fig. 3. — Same data as for figure 2, but with model spectra for pure thermal and nonthermal 
bremsstrahlung corresponding to an electron Alfven velocities of 56, 68, and 80 atomic units, for 
fully ionized helium. 
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